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Molecular chirality leads to a wonderful variety of equilibrium structures, 
from the simple cholesteric phase to the twist-grain-boundary phases, and 
it is responsible for interesting and technologically important materials 
like ferroelectric liquid crystals. This paper will review some recent 
advances in our understanding of the connection between the chiral 
geometry of individual molecules and the important phenomenological 
parameters that determine macroscopic chiral structure. It will then 
consider chiral structure in columnar systems and propose a new 
equilibrium phase consisting of a regular lattice of twisted ropes. 
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1 Introduction 

Chirality leads to a marvelous variety of liquid crystalline phases, including 
the cholesteric, blue, TGB, and SmC* phases, with characteristic length 
scales in the "mesoscopic" range from a fraction of a micron to tens of mi- 
crons or more. In spite of its importance, remarkably little is known about 
the connection between chirality at the molecular level and the macroscopic 
chiral structure of liquid crystalline phases. In particular, it is not possible 
at this time to predict the pitch P (of order a micron) of a cholesteric phase 
from the structure of its constituent molecules. This talk will address some 
aspects of the molecular origins of chirality. It will also speculate about a 
possible new phase in chiral columnar or polymeric systems. 
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2 What is Chirality? 

A molecule is chiral if it cannot be brought into coincidence with its mir- 
ror image]]]]. Examples of chiral molecules and the achiral configurations 
from which they are derived are shown in Fig. [[]. A chiral molecule must 
have a three dimensional structure. A linear, "one-dimensional," or a flat, 
"planar," molecule cannot be chiral. In particular a chiral molecule cannot 
be uniaxial. Nevertheless, the local structure of liquid crystals phases (in- 
cluding the cholesteric and other chiral phases) is nearly uniaxial. Indeed 
the phenomenology of most chiral phases can be explained in terms of the 
chiral Frank free-energy density, 

/= iK 1 (V-n) 2 + |ir 2 [n-(Vxn)] 2 + iK 3 [nx(Vxn)] 2 -/in-(Vxn), (1) 

expressed in terms of the Frank director n(x) specifying the local direction 
of uniaxial molecular alignment at the space point x. This free energy 
does not distinguish directly between truly uniaxial molecules and biaxial 
molecules that spin about some local axis so that their average configu- 
rations are uniaxial. Chirality in this free energy is reflected in the term 
—hn ■ (V x n). The parameter h is a chiral or pseudoscalar field that is 
nonzero only if the constituent molecules are chiral. It is a phenomeno- 
logical parameter that may vary with temperature (or other external field 
such as pressure) and may even change sign. Note that the Frank free en- 
ergy makes no explicit reference to the fact that a chiral molecule cannot 
be uniaxial, h determines the length scale of equilibrium chiral structures. 
For example, the pitch wavenumber ko = 2tc/P of a cholesteric is simply 
— h/K 2 . Thus the chiral parameter h can be identified with K 2 k . It has 
units of energy/ (length) 2 , and one would naively expect it to have a mag- 
nitude of order T NI /a 2 , where T NI is the isotropic-to-nematic transition 
temperature and a ~ 30A is a molecular length. The elastic constant K 2 
does conform to naive expectations from dimensional analysis and is of or- 
der T NI /a. Thus h = (2itK 2 /P) « 2nT NI /(aP) is or order [T NI /a 2 ](a/P), 
which is a/P ~ 10 -3 times smaller than naive dimensional analysis would 
predict. 

The Frank free energy of Eq. (|I]) can be used to treat the cholesteric 
phase, chiral smectic phases (when a smectic free energy is added) and 
even low-chirality blue phases. To describe transitions from the isotropic 
phase to the cholesteric or blue phases, the Landau-de-Gennes free energy 
expressed in terms of the symmetric-traceless nematic order parameter Q l i 
is more useful. When there are chiral molecules, this free energy is 

f Q = |rTrQ 2 - wTiQ 3 + uTiQ 4 

+lC 1 d i Q> k d i Q> k + ICzdiCtdjCy* + HeWCfdjQ". (2) 
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Figure 1: Examples of chiral structures created from achiral ones: (a) left, 
achiral tetrahedron with four equal masses at its vertices; right, a chiral 
tetrahedron with four unequal masses at its vertices, (b) left, an achiral 
planar "H"; right, a chiral "twisted H". (c) left, an achiral propeller with 
all blades perpendicular to the hexagonal core; right, a chiral propeller with 
all blades rotated away from the normal to the hexagonal plane, (d) left, 
an achiral planar sheet; right, a helix formed by twisting a sheet about an 
cylinder, (e) left and middle, cylinder with achiral deorations; right, chiral 
cylinder with helical decorations. 
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The chiral term in this expression is the one proportional to the Levi-Civita 
symbol e* jfc . Indeed, the usual mathematical signal of chiral symmetry 
breaking in a free energy is a term linear in e* jfc . It is present in the Frank 
free energy in the term n • (V x n) = e i: ' k n l djn k . In the nematic phase the 
order parameter is unixial: = S{n l n^ — |<5 i -'), where S is the Maier- 
Saupe order parameter, and Eq. (|J) reduces to Eq. ([I]) with h ~ HS 2 . 
Thus, like the Frank free energy, the Landau-de-Gennes free energy makes 
no explicitly reference to biaxiality. 

A theoretical goal should be the calculation of the parameter h (or 
H) from molecular parameters and inter-molecular interactions. A the- 
ory should (1) show why h is smaller than naive arguments would pre- 
dict, (2) provide at least semi-quantitative guidance as to how molecular 
architecture affects h, and (3) elucidate the transition from necessarily 
non-uniaxial chiral molecules to an essentially uniaxial macroscopic free 
energy. Since h is zero for achiral molecules, it is natural to expect that 
h will be proportional to some parameter measuring the "degree of chiral- 
ity" of constituent molecules. Section III will introduce various measures 
of molecular chirality and discuss why it is impossible to define a single 
measure of chirality. h should depend on interactions between molecules 
as well as on chiral strength. Section IV will present a calculation of h for 
simple model molecules composed of atoms interacting via classical cen- 
tral force potentials 0. This calculation will highlight the role of biaxial 
correlations. It ignores, however, chiral quantum mechanical dispersion 
forces 0. Finally, Sec. V will discuss chirality in columnar systems and 
speculate on the possibility of a new chiral phase consisting of a lattice of 
twisted ropes. 



3 Measures of Molecular Chirality 

To keep our discussion as simple as possible, we will consider only molecules 
composed of neutral atoms. All information about the symmetry of such 
a molecule can be constructed from its mass density relative to its center 
of mass: 

p(x) = ^m a 5(x - r a ), (3) 

where m a is the mass of atom a whose position vector relative to the 
molecular center of mass is r a . A molecule is chiral if there exists no 
mirror operation M under which p(x) is invariant, i.e., for all M, p(x) ^ 
p(Mx). Mass moment tensors contain information about the symmetry 
and mass distribution of a molecule, and they can be used to construct 
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chiral measures. The first moment specifies the position of the center of 
mass, which is uninteresting since we can always place the center of mass 
at the origin of our coordinates. Second- and third-rank moments are, 
however, interesting, and we define 



a 



(4) 



and 



D kim = J2m a [r k a r k y a -l(r a ) 2 (r k J lm + r l J km + r™5 kl 

a 

D k n lm = E^(^r[^i-i(^) 2 (^ m +^ fem +coj- (5) 

where (r a ) n is (r Q • r a ) n l 2 . These tensors are symmetric and traceless, 
and they are invariant under the mirror operation that interchanges any 
two axes, say x and y. Therefore, they cannot by themselves encode any 
information about molecular chirality. 

To construct a quantity that is sensitive to chirality, we need com- 
binations involving three distinct tensors in combination with the anti- 
symmetric Levi-Civita symbol e^ k . One possible chiral parameter is 
e tjkjjiim(jji^(j2^km ^ w h ere ^he Einstein convention on repeated indices is 
understood, and where (C 2 ) km is the km component of the tensor C 2 . 
More useful chiral measures can be constructed by dividing C* J into a uni- 
axial and a biaxial part. The tensor C l i has five independent components 
that can be parametrized by two eigenvalues ipQ and tpB and an ortho nor- 
mal triad {ei,e2,e3} specifying the directions of the principal axes of the 
molecule. Thus, we can write 



where 
and 



C ij = Qij + B 



= Me\e{-elei)=tp B B 



(6) 

(7) 
(8) 



The last equation defines the reduced biaxial tensor B lj that depends only 
on the principal axis vectors ei and e2. For a uniaxial molecule, the pa- 
rameter ips is zero, and we will refer to Q 1 ^ as the uniaxial part of C 1 - 7 
and B l i as the biaxial part (even though the equilibrium average ) can 
develop a biaxial part). The third rank tensor has seven independent 
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components, which can be represented in the basis {ei, 62,63} defined by 
the principal axes of C ij . Of these, only one component is of interest to 
the present discussion: 

D klm = ^ c {e\e\e™ + 5 perm) = ip c D klm \ (9) 

This equation defines the reduced third-rank tensor £) y . We can now 
introduce the chiral strength 

tj) = e ijk eie l 3 B jm D klm = 2ip B ip D = D xyz (B xx - B yy ). (10) 

This parameter clearly changes sign under any mirror operation. It must, 
therefore, be zero for any achiral molecule. It is a continuous function 
that has opposite signs for structures of opposite chirality. Consider, for 
example, the following mirror operations: (1) z — > —z, D xyz changes signs 
but B xx — B yy does not; (2) x <-> y, (B xx — B yy ) changes sign but D xyz 
does not. More generally, we could introduce a class of chiral strength 
parameters defined in terms of D klm and the biaxial part, B^ , of the tensor 

W = ^^t' m (11) 
The parameter -?/>oo is equal to ip. For all other p and n, ip pn is distinct from 

It is instructive to calculate ijj for some of the chiral molecules shown in 
Fig. |l| The "twisted H" structure is the simplest example of a geometrically 
chiral object. It has the pedagogical virtue that its geometric properties 
depend in a simple continuous way on the angle 7 between arms. When 
7 = 0or7 = 7r/2, the H is achiral. Thus any chiral measure must be 
zero at these points, and our expectation is that chiral measures will go 
continuously to zero at these points. Let r be the radius of the arms of 
the "H" and d be half its height. Then ip Q = 2(2d 2 - r 2 ), ip B = 2r 2 cos 7, 
ipD = 2c?r 2 sin7, and 

ip = 4dr 4 sin27. (12) 

As required, ip — > as 7 — > 0, 7r/2. Another interesting example is the helix 
or spiral structure in Fig. |T](d). Parametrizing the helix by its position as a 
function of normalized arclength s: R(s) = [r cos(2ims),r sin(2ims), Ls], 
where n is an integer, we find 

3r 4 L 



(2nn) 



24 



(7m) 2 



(13) 



Thus, ip decreases as the helix becomes more tightly wound (n increases) - 
a reasonable result because at large distances a tightly wound spiral looks 
more like a uniform cylinder than does a not so tightly wound one. 
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There are many, in fact infinitely many, chiral strength parameters that 
provide a measure of the degree of chirality of a molecule. The set ip pn de- 
fined above is just one of infinitely many sets we can define. At first, it 
may seem strange that the degree of chirality may be characterized in so 
many different ways. After all, chirality is simply the absence of a mirror 
symmetry, which involves a discrete operation. We might, therefore, have 
expected a single "chiral order parameter" analagous to the order parame- 
ter for an Ising model, which also describes the absence of a discrete sym- 
metry. Chirality is, however, a function of the positions of all of the atoms 
in a molecule, which cannot be characterized by a single variable. Thus, 
chirality is perhaps more like the absence of spherical symmetry. There an 
infinite number of parameters, the spherical harmonics for example, that 
chracterize deviations from perfect sphericity. 

Some measures of chirality are more useful than others, As we will see, 
the large- distance potential between two generic chiral molecules is directly 
proportional to the strength ip introduced in Eq. (|iTj|). Thus, if) provides 
a good measure of the strength of chiral interactions and enters into the 
calculations of h. There are molecules, such as the chiral tetrahedron of 
Fig. |I|a, for which if) is zero. For such molecules, if) is clearly not very 
useful. Other measures are, however, nonzero and provide a measure of 
the strength of chiral interactions. 

4 Chiral Pathways between Enantiomers 

The mirror image of a chiral molecule is called its enantiomer. Any measure 
of chirality will change sign under a mirror operation, and as a result, chiral 
enantiomers will have chiral strengths of equal magnitude but opposite 
sign. This property implies that any series of continuous distortions of a 
molecule that takes it from a chiral configuration to that of its enantiomer 
must necessarily pass through a point where any given chiral measure is 
zero. It is tempting to conclude that this point corresponds to an achiral 
configuration and that any pathway between chiral enantiomers must pass 
through an achiral state. This conclusion is false. There are continuous 
pathways between chiral enantiomers all of whose states are chiral (for 
further discussion of this and the concept of a "topological rubber glove," 
see . A molecule is achiral if and only if any and all chiral measures are 
zero. For a particular chiral configuration, a given chiral measure may be 
zero, but at least one chiral measure be must be nonzero. 

To illustrate this point, let us consider the pathways depicted in Fig. 
between reflected images of a twisted H. A is a twisted H with an angle 
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A F E 




BCD 

Figure 2: Pathways between chiral enantiomers consisting of 4 atoms at 
the vertices of an H. A and E are chiral enantiomers (mirror images). 
The 4 masses and their positions can be changed continuously to provide 
continuous pathways between the two enantiomers. Path AFE passes 
through the achiral planar configuration F. Path ABCDE passes through 
chiral configurations only. Configuration C is the chiral tetrahedron with 
unequal masses. 
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7 = 7r/6 between top and bottom arms. E is its chiral enantiomer with 
angle —7r/6 between arms. F is the the achiral planar structure. B through 
D are all chiral configurations with different masses at the vertices of the 
H. C is the chiral tetrahedron with four unequal masses. The path AFE 
between chiral enantiomers passes through the achiral configuration F for 
which all chiral measures are zero. On the other hand, all configurations in 
the path ABCDE are chiral even though A and E are chiral enantiomers 
for which any chiral measure ip must have opposite signs: ip{E) = —ijj(A). 
Thus, any continuous chiral measure must pass through zero at some point 
on the path ABCDE, but different chiral measures do not have to pass 
through zero at the same point. To illustrate these facts, we consider 
a specific parametrization of the closed path ABCDEFA. Let the four 
masses in the twisted H be 

{mi, m 2 , m 3 , m 4 } = m[{l, 1, 1, 1} + /x{0.31, -0.31, 0.23, -0.23}]. (14) 

The configurations of the twisted H are thus specified by the three param- 
eters /i, d/r, the ratio of the semiheight to the arm radius, and 7, the angle 
between the arms of the H. The following properties of a chiral measure ip 
are apparent: 

1. fi = 0: Equal masses 

(a) 7 = 0, 7 = 7r/2: achiral =^■0 = 

(b) d = rsin(7/2) or d — rcos(7/2): achiral =>- ip = 

(c) = ~4>{~l) = — "^(tt — 7) = 2c/r 4 sin27. 

2. /i 7^ 0: Unequal masses 

(a) 7 = 0: planar achiral molecule =>■ ip = 

(b) 7 7^ 0: chiral molecule =>- ip 7^ 

(c) = ^(7,-j") 

(d) d = r I chiral tetrahedron with unequal masses. 

The points in Fig. ^ correspond to the following values of [7, fi, d/r]: 
A = [vr/6,0,1], B = [tt/3,3/4,(V2 + 3)/(4V2)], C = [0,1,1/V§], D = 
[271-/3,3/4, (V2 + 3)/(4V2)], E = [5vr/6,0, 1] = [-tt/6,0,1]. Figure |a 
shows the normalized chiral strengths \l/ nm = z/w(7> £*> d/r)/ij){A) for 
(n, m) = (0, 0) and (n, m) = (0, 8) as a function of 7 for the path EFA with 
fi = 0, d/r = 1. Only one curve appears in the figure because the two curves 
for ^00 an d ^08 are identical. As required, the curve passes through zero at 
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the achiral point 7 = 0. Figures |3|a and [5]b show the normalized strengths 
* = \[r 00 and \& 8 for the path ABCDE with 7r/6 < 7 < 57r/6 with 
/i = 1 - (9/4) [(2 7 /tt) - 1] 2 and d/r = (1 V2) + (9/4) [1 - (1/V2)] [(2 7 /tt) - 1] 2 . 
The two curves are clearly different, but they appear to go to zero at 
7 = 7r/2. Figure ^|d shows a blowup of the region around 7 = 7r/2, showing 
that ^ and \l/os pass through zero at different values of 7. Note that \I/ 
passes through zero at 7 = 7r/2 indicating that this parameter does not 
provide a measure of the chirality of the tetrahedron with unequal masses. 
\l/o8 is nonzero but small at 7 = n/2. 

The path ABCDE was chosen to pass through the chiral tetrahedron 
at point C. A tetrahedrally coordinated carbon atom bonded to four dif- 
ferenct atoms or chemical groups is a chiral center. It is used to identify 
chriality in molecules. Interestingly, the simple chiral tetrahedron is not 
very chiral according to the measures ip pn . All measures ip 00 to tp 05 are 
zero for this molecule. This is a partial explanation of why the difference 
between \1/ and \l/os is so small near 7 = n/2. Other pathways between A 
and E would display more pronounced differences. 

5 Calculation of h 

Two distinct microscopic mechanisms contribute to a non- vanishing macro- 
scopic chiral parameter h. The first is quantum mechanical in origin and 
has its clearest manifestation in the generalization of the familiar Van der 
Waals dispersion potential to chiral molecules The electrostatic poten- 
tial between molecules is expanded to include dipole-quadrupole as well 
as dipole-diple interactions. The expectation value of the dipole-moment- 
quadrupole- moment product (p t Q : ' k ) is nonzero for a chiral molecules and 
is proportional to (e^e^e* + perm) for a molecule spun about the direction 
e 3 . The resultant potential has the form, 

U = V(R)e 3 -e' 3 e 3 -(Rxe' 3 ), (15) 

where R is the spatial separation of the centers of mass of the two molecules 
and V(R) ~ R~ 7 . This potential leads directly to the Landau-de-Gennes 
chiral interaction Q &*Q in the isotropic phase [Eq. (0)] and to the 
Frank chiral term in the nematic phase with h ~ S 2 . The calculation of 
the strength of the potential V(R) for real liquid crystal molecules is com- 
plex. The second mechanism can be described in purely classical terms. It 
arises from central force potentials, including hard-core potentials, between 
points (or atoms) on molecules with chiral shapes. This is the mechanism 
popularized by Straley'sf| image of two interlocking screws that twist rel- 
ative to each other. Here we will consider only the second mechanism. 
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Figure 3: Plots of reduced chiral strength parameters for different path- 
ways between chiral enantiomers shown in Fig. |2|. (a) $ = ipoo/ipoo(A) = 
ipos / i[)08{A) for the path EFA passing through the achiral configuration 
F. (b) = W<M^) for the path ABODE, (c) ^ 08 = WV>os(^) for 
the path ABODE, (d) \I/ and \l/o8 i n the vicinity of 7 = n/2. They pass 
through zero at different points. 
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We will show that h is nonzero only if there are biaxial correlations be- 
tween molecules and as a consequence is zero in any mean-field theory with 
uniaxial order parameters. 

We consider each molecule A to be composed of atoms a whose positions 
within the molecule are rigidly fixed. Atoms on different molecules interact 
via a central potential V(R). Let Ha be the position of the center of 
mass of molecule A, be the position of atom a relative to R^, and 
R^a = Ra + r Aa be the position of atom a in molecule A. The total 
intermolecular potential energy is thus 

U = \ E V(\R Aa -R B „\). (16) 

AajiB0 

The chiral parameter h is the derivative of the free energy density with 
respect to the cholesteric wave number k = 2n/P evaluated in the aligned 
state with spatially uniform director n and k — 0: 



where Q is the system volume and 




= -^T.Ts„ (17) 

fc=0 ^" BA 



Tba = (JZ^ kR Hd k V{R + vp-v a )) (18) 

\ (3a I 

is the "projected torque", R^ • tab, where tab is the torque exerted on 
A by B. Here r = R# — R^ and is the component of r perpendicular 
to the uniform director n. This expression for h is quite general. It can 
be used in Monte Carlo and molecular dynamics simulations. It is valid 
both classically and quantum mechanically. (It will reproduce the disper- 
sion calculations described above if electron as well as atomic coordinates 
are included). Here we restrict our attention to classical central forces. 
The fundamental question is, how do nonchiral central forces give rise to 
effective chiral forces? To answer this question, we consider the multipole 
expansion of U . Schematically, we have 

J2V(R A a-RB(s) - l[d^d k d l V}C%C k B l (19) 

+^[d l d j d k d l d m V}(C i iDf n + C^Df"). 

The second term on the right hand side of this equation is the chiral term. 
It can be expressed as 

V c = y> D il) B [&d i &&d n V\(I?iD'£ m + B%D k A lm ), (20) 
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where the reduced tensors B l i and D klm are defined in Eqs. (|8|) and (|9|). 
The triad of vectors {ei, 62,63} forms a right-handed coordinate system. 
Since this triad is determined by the second-moment tensor, it can be cho- 
sen to be invariant under a chiral operation M. Thus, the reduced tensors 
B^ and D klm do not change sign under a chiral operation. The chiral pa- 
rameter ip = 2ipB^D and the potential d 5 V both change sign under a chiral 
operation on all atoms. Thus the potential Vc is invariant under a chiral 
operation involving all atoms. If, however, we perform a chiral operation 
about the center of mass of each molecule (i.e., convert each molecule to 
its enantiomer) leaving the centers-of-mass positions fixed, ip changes sign 
but d 5 V remains unchanged. Thus, the potential Vc changes sign upon 
reversing the chirality of constituent molecules. As we noted earlier, any 
term linear in the Levi-Civita symbol e u ' fc is chiral. The potential Vc is lin- 
ear in e l i k because the triad {ex, e 2 , e 3 } is right handed, and ex = (e 2 x e 3 ), 
^2 = (^3 x ex), and e 3 = (ex x e 2 ). Thus 

D ijk = e ilm eie[B km + t klm ^e[B jm + e jlm e k e l 3 B im . (21) 

We now consider the evaluation of h from Eq. (|T7|) and the multipole 
expansion of Eq. fl2"0|). To simplify our calculation, we restrict all molecular 
axes 63,4 to be rigidly aligned along a common axis no- This is equivalent 
to forcing the Maier-Saupe order parameter S to be unity. In this case, 
the potential Vc involves only the product of the biaxial moment B^ on 
different molecules, and upon averaging, the torque becomes 

T BA = ^(K(R)T B (R)), (22) 

where 

T B (R) = (B%B l l) = 2(cos2(0 B - cj> A )) (23) 
is the biaxial correlation function and 

K(R) = R\ {g^ + R\g^ + \Rl[g^ + R*g®}} , (24) 

where g(R 2 /2) = V{R), and g {n) {x) = d n g(x)/dx n . If molecules rotate 
independently, the biaxial correlation function r^(R) is zero. Mean- field 
theory seeks the best single particle density matrix and thereby ignores 
correlations between different particle. Thus, in mean-field theory, is 
zero, and the chiral potential h is zero unless there is long-range biaxial 
order with a nonvanishing value of (cos 20^). The latter observation was 
made some time ago|| 0. Nevertheless, there have since been mean-field 
calculations claiming to produce a nonvanishing h in a uniaxial system§|. 
In uniaxial systems, Tb(R) oc e~ R ^ where £ is the biaxial correlation, which 
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should be order a molecular spacing. We can now provide a partial answer 
to the question of why h is so much smaller than naive scaling arguments 
would predict. First, the chiral strength ip may be quite small, particularly 
for tightly would helical molecules such as DNA or tobacco mosaic viruses. 
Second, the biaxial correlation length could be small (because, for example 
the molecules are only weakly biaxial), allowing only very near neighbors 
to contribute to the sum over AB in Eq. (|l7|) for h. 



6 Chiral Columnar Phases 

We now have a fairly complete catalog and understanding of equilibrium 
phases of rod-like chiral mesogens. Our understanding of chiral phases 
of disc-like chiral mesogens is less complete. Such mesogens have been 
synthesized P, [T0|1 . They produce cholesteric and columnar mesophases, the 
latter of which exhibit unusual spiral textures 0. Columnar phases with 
spontaneous chiral symmetry breaking in which propeller-like molecules ro- 



tate along columns have been observed [IT]. In addition, long, semi-flexible 
polymers like DNA whose natural low-temperature phase has columnar 
symmetry, exhibit cholesteric, columnar, and blue phasesjl^], and even 



what appears to be a hexatic phase|13|]. There have been theoretical pre 



dictions of a number of phases including (1) a moire phaseJT^] in which the 
hexagonal columnar lattice rotates in discrete jumps about an axis parallel 
to the columns across twist grain boundaries consisting of a honeycomb 
lattice of screw dislocations, (2) a tilt-grain boundary phase [|nj in which 
columns rotate in discrete jumps bout an axis perpendicular to the columns 
across twist-grain boundaries consisting of a parallel array of screw disloca- 
tions, and (3) a soliton phase [15| in which molecules all rotate in the same 
sense in their columns without disrupting the regular hexagonal columnar 
lattice. Many other states are possible. Here we will consider one of them: 
a hexagonal lattice of twisted ropes. This phase would be constructed as 
follows: First, take cylindrical sections of radius R of a columnar lattice 
and twist them about an axis parallel to the columns as shown in Fig. f|. 
This operation will cost strain energy but will gain twist energy. Next put 
these columns on a regular hexagonal lattice and deform the cylinders to a 
hexagonal shape to produce a dense packing as shown in Fig. [5]. The inter- 
face between adjacent hexagonal cylinders is a twist grain boundary and 
is favored if chiral interactions are sufficiently strong. On the other hand, 
there is a strain energy cost associated with distorting the cylinders and 
an energy cost associated with the vertices of the hexagons. The radius of 
the original cylinders (and thus the lattice spacing) and the degree of twist 
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Figure 4: Twisting a cylinder of a columnar lattice, (a) The untwisted 
cylinder of radius R with equilibrium spacing between columns d. (b) The 
twisted cylinder, (c) The outer surface of the twisted cylinder unrolled. 
The columns make an angle 9 = tan -1 k R with respect to the cylindrical 
axis. The normal spacing between between layers is now dj cos 9 « d(l + 
9 2 /2). 



within a cylinder will adjust to minimize the free energy. This lattice of 
twisted ropes should be competitive with the moire, tilt-grain-boundary, 
and soliton phases. 

To be more quantitative, we now investigate in more detail the various 
contributions to the energy of this lattice. The elastic energy density of a 
columnar phase with columns aligned along the z axis is 



where 



/ei = |A«L + fiu af3 u a/3 + \Kz{d x u) 2 + §K 6 (W 6 ) 2 , (25) 



u a /3 = UdaUpdpUa + d^u a d^Ufi - d z u a d z up), a, (3 = x,y (26) 



is the nonlinear strain for a columnar lattice in mixed Euler-Lagrangian 
coordinates. The first three terms in this nonlinear strain are identical to 
the nonlinear Lagrangian strain (with a relative plus sign between the linear 
and nonlinear terms) of a two-dimensional solid in which the coordinates 
x and y refer to positions in a reference unstrained sample. The third 
term reflects the three-dimensional nature of the columnar phase. There 
is a minus sign between it and the preceding terms because z measures 
a coordinate in space and not a coordinate in a reference solids; it is an 
Eulerian rather than a Lagrangian variable jl6| . The variable 9q = (d x u y — 



d y u x )/2 in Eq. ( [25] ) specifies the direction of hexatic order of the lattice. In 
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Figure 5: Close-packed twisted cylinders deformed to spacefilling hexagonal 
twisted cylinders. 

a chiral system, there are additional chiral terms in the free energy density: 

/ch = -led z 9 6 - 7 „n • (V x n), (27) 

where n is the nematic director, which is parallel to the local columnar axis. 
The first term in f c ^ favors twisting of the hexagonal columnar lattice along 
the columns, i.e., it favors a twisted rope or the moire configuration. The 
second term favors twisting of the columns along directions perpendicular 
to the local columnar axes, i.e., it favors a tilt-grain-boundary structure. 

We can use / e i and / c h to calculate the preferred configuration of a 
cylinder of a chiral columnar lattice of radius R in which dislocations and 
other defects are not allowed. The elastic energy of a lattice in which 
# 6 = k z is uniformly twisted with twist wavenumber k is controlled by 
the nonlinear term d z u a d z up of the strain. Indeed the two-dimensional 
Lagrangian part of the strain arising from such a uniform twist is zero! 
There are nonlinear contributions to the strain, however, because the dis- 
tance between neighboring columns increases as the columns are twisted 
as depicted in Fig. |j. Since twisting produces no linear strain, the energy 
required to twist a columnar lattice is lower than that required to twist 
a three-dimensional elastic solid. The elastic energy per unit length of a 
twisted columnar lattice can be calculated. For small k R, it is 

-p = JL(k R) 4 7iR 2 + lK 3 kl(k R) 2 7iR 2 + \KsklvR 2 . (28) 

LJ 

where ~p is an effective modulus that is strictly proportional to fi in an 
incompressible system where the bulk coefliceint A is infinite. This energy 
should be compared to the result, fik^R 4 , for a three-dimensional solid. 
The chiral energy per unit length is 

^ = -bok + ln k (k R) 2 }nR 2 . (29) 
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Figure 6: Creation of a twist grain boundary at the interface between two 
cells. All cells twist in the same direction. The columns in back (-) and in 
from (+) of a given interface twist in opposite directions. The interface is 
thus a twist grain boundary consisting of screw dislocations perpendicular 
to the z axis separated by a distance Id = d/2 sin 6, where 29 is the change 
in angle of the columnar lattice across the boundary. 

Equations ( P%D and (^) can be used to determine the twist wavenumber 
ko for a cylindrical rod of radius R. 

Our goal is to determine the lowest energy lattice of rods for given values 
of the coupling constants. We begin by distorting our cylindrical rods into 
hexagonal rods so that they can pack closely in an hexagonal array (Fig. |5|). 
Each unit cell of the lattice is now a hexagonal cylinder of twisted columns. 
This distortion will cost an energy for each rod proportional to its shear 
modulus times its volume irR 2 L: -Edist. = afinR 2 L, where a is a constant 
of order unity. The columns at the outer edge of the hexagonal unit cells 
make an angle of 9 = tan -1 k$R « koR relative to the vertical z axis. 
Columns on opposite sides of a boundary between adjacent hexagonal unit 
cells, however, make opposite angles relative to this axis, so the boundary is 
in fact a twist grain boundary across which the columnar lattice undergoes 
a change in angle of 26 as shown in Fig. |6|. There is a positive energy 
cost to a grain boundary coming from the cost of creating the dislocations 
in the boundary. There is also an energy gain in chiral systems arising 
from the — 7 n n ■ (V x n) term in f c ^. Consider a boundary in the xz 
plane with dimensions L x x L z . Since J dyn • V x n = L X L Z 28, where 
the integral is across the grain boundary, there is a chiral energy gain of 
—•y n L x L z 29 ~ —•y n L x L z 2koR. The dislocation energy is -Edisc = NdL x e, 
where e is the energy per unit length of a dislocation, and Nj = L z /ld is 
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the number of dislocations with l d the distance between dislocations. There 
is a simple geometric relation between Id and 9: 9 = sirT 1 d/2ld ~ d/2ld, 
where d is the spacing between columns in the columnar lattice. Thus the 
energy of a grain boundary for small k$R is Eqb = L x L z {ed~ x — 2 r y n )2koR. 
The length L x is proportional to R: L x = anR/3, where a is a constant of 
order unity. The first term in / c h contributes an energy —TcR 2 L z ^ekQ per 
column. Finally, there is an energy per unit length r/2 for each vertex of 
the hexagonal lattice where grain boundaries meet. Combining all of these 
energies and keeping only the dominant terms, we find 

/ = f = HkoR) A + ^ + an, (30) 

where 7 = 70 + (4a/3)(7„ — erf -1 ). Minimization of this energy over k for 
7 > yields 




^0 = h=W ~ R~" ■ (3i) 



Thus k R ~ R x / 3 is much less than one when 7 is small and R is large. 
At this optimal value for k , 

^-H^r vi/3 +w < 32 > 

Minimization over R yields 

*-e&£ ,33, 

and 

1 7 4 

/ = -7T=^+«/i. (34) 

Thus, for sufficiently strong chiral coupling 7, this free energy will be lower 
than that of a uniform columnar phase with no twist, and it is possible 
that there is some range of parameters for which this phase has a lower 
energy than either the moire or the soliton phase. 
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